NASA Technical Memorandum 89494 


Results of the Integration of 
a Transonic Full-Potential 
Analysis Program with a Free- 
Wake Lifting-Line Program for 
Hovering Rotors 

Song-Young Chung, Ames Research Center, Moffett Field, California 

<NASA-TH-89494) BESOIIS C t IBB INSBGBATIOH N87-22621 

Cf A IE A li SC NIC fUIl-PCTENTIAL ANALYSIS 
EEOGEAH HUH A EBEE-HAKE IIEIING-LIHE 

fBCGBAfl I OB HCVEBING EOTOHS (NASA) 22 p <3.3 Oaclas 
Avail: Nils BC A02/BE A01 CSCL 01A *1/02 0076844 

May 1987 


NASA 

National Aeronautics and 
Space Administration 


Ames Research Center 

Moffett Field. California 94035 


NASA Technical Memorandum 89494 


Results of the Integration of 
a Transonic Full-Potential 
Analysis Program with a Free- 
Wake Lifting-Line Program for 
Hovering Rotors 

Song-Young Chung 


May 1987 


rw\s/\ 

National Aeronautics and 
Space Administration 




SUMMARY 


This report discusses the hovering performance predictions of the TFAR1 and 
OPLIN codes and the experimental data discussed in reference 6. The TFAR1 program 
solves the full -potential equation in a rotor-fixed coordinate system by use of the 
line relaxation method. The OPLIN program calculates the positions of wake vortices 
and rotor performance using the influence-coefficient-and-lifting-line method. The 
two programs are combined by adding the induced velocities from the OPLIN code to 
the near flow-field of the rotor from the TFAR1 code. Results of this investigation 
show that the TFAR1 program converges better with the downwash-coupling method than 
the twist correction method to include the wake downwash. 
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Recently, several finite-difference methods using small-disturbance 
(refs. 1,2), full-potential (refs. 3,4), or Euler equations (ref. 5) have been used 
for the analysis of helicopter rotor flow. Wake effects are computed by either 
momentum theory or lifting surface methods. Experimental and analytical results for 
hovering rotors are discussed in reference 6 where the tip vortices were found to 
reach their maximum strength very quickly, but not in the core size or radial and 
axial positions. The best agreement between experiment and analytical results was 
obtained when the tip vortex contraction was reduced by 0.025 R (radius) from the 
experimental data. 

After the effective angle of attack has been determined from momentum theory 
and blade flapping motion, the pressure distribution on a helicopter rotor was 
calculated in reference 1 from the small-disturbance potential equation. In refer- 
ence 2, the effects of transonic flow on the advancing blade are calculated. by using 
an unsteady, transonic small-perturbation method with a Drees downwash model. In 
reference 3, the full-potential equation is derived using body-conforming coordi- 
nates fixed on the rotor, and the pressure distributions computed on the advancing 
blade are presented for nonlifting helicopter rotors with sweep and taper. The 
TFAR1 program used in the present report was developed in reference 3 where the one- 
dimensional ( 1— D) and 2-D wake models were introduced. Lifting calculation for the 
advancing blade were presented in reference 4 where the effective angle of attack 
was calculated from the wake portion of the program CAMRAD (ref. 7) at each spanwise 
station. 

In reference 5 a semi-implicit approximate factorization is used to factorize 
the unsteady Euler equation where the spanwise derivative is lagged behind the other 


1 



two derivatives by one time-step. This semi-implicit Euler equation is solved for 
wing and hovering rotors. Various finite-difference and lifting-surface methods, 
using acceleration or velocity potential, are reviewed with wake modelling in refer- 
ence 8 for rotary wing analysis. The capability of the Euler code to capture wake 
vortices is discussed. 

In this paper, the free-wake lifting-line program OPLIN, as developed in refer- 
ence 9 and 10-14, is coupled with the program, TFAR1, for the analysis of hovering 
rotors with transonic tip speeds. 

The author wishes to express thanks to Dr. I-Chung Chang for allowing me to 
couple his code, TFAR1, with my code, OPLIN. 


FULL - POTENT I AL FORMULATION 


The velocity components in a rotating coordinate reference frame are decomposed 
into the freestream component, perturbation potential, and the additional downwash 
caused by the wake vortices which were not captured by the transonic full -potential 
equation. The velocity components are 
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where q is the total velocity, v is the velocity caused by the rotor rotation, 

4 > is the perturbation potential, w is the downwash calculated everywhere by the 
free-wake lifting-line code in the computational domain, and Q is the rotor rota- 
tional velocity. 

The circulation r at each span station is defined by the potential jump at 
the trailing edge. 


r = q • dr = V$ • dr = d$ = * u - (7) 

The potential jump is the wake vortex strength in the 3-D flow. The conservation of 
circulation for inviscid flow and the force-free condition on the wake are satisfied 
when the potential jump is convected with the local velocity. The Kutta condition 
defines the potential jump at the trailing edge to be the magnitude of the 
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circulation of the lifting body. This jump is interpolated into the numerical grid 
line as follows. 

*1u = *3u + = *31 + A * (8) 

*11 = *3u • t*l = *3u " A * (9) 

Here A represents the jump in the quantity. Also, the flow-tangency condition on 
the body is applied according to the following numerical boundary conditions. 

q • n = + (V + w) =0 (10) 

♦n n 

^ + 2An(V + w) n = 0 (11) 

♦l = $3 + 2An(V + w) n = 0 (12) 

FREE WAKE FORMULATION 


The free-wake lifting-line theory developed in reference 8 divides the wake of 
a hovering rotor into a near wake attached to the blade trailing edge, an interme- 
diate wake of several inner and tip ring vortices, and the far wake of semi-infinite 
vortex cylinders (fig. 1). The downwash induced by the entire wake on the rotor 
blade and in the wake vortices is computed as influence coefficients which are a 
function of wake geometry. 


On the position i of the quarter-chord line of the blade, or in the interme- 
diate wake, the induced velocities are (w ^ + v i) , such that 
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where y is the strength of the bound circulation, and 
y- and z- direction influence coefficients at position 
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A . . are the 
21, J 


yi.J 

i caused by all of the 
trailing vortices of bound circulation j of unit strength. The influence coeffi- 
cients are computed by using the Biot-Savart law from the wake geometry of the 
v th iteration. The bound circulation of the (v + 1)th iteration is computed by 
applying the Kutta-Joukowski law at each blade section. The wake geometry of the 
(i + 1)th position is computed by integrating the downwash on the wake positions 
at i and i - 1 as follows. 


z i = z t-1 * °- 5 * (w i-1 

+ W.)* 27T/b 
1 

(14) 

r i = r i-l * °- 5 * (v i-1 

+ v.)* 2ir/b 

( 15) 


where z and r are the axial and radial positions of the intermediate and far-wake 
vortices, w and v are the corresponding induced velocities, and b is the number 
of blades. The r and z for i = 0 are the radial and axial positions of the near 
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wake considered to be rolled up instantaneously on the plane of rotation according 
to the law of conservation of linear momentum as follows. 

(r k - Vo * J * , P j 4r j (16) 

The bound circulation is calculated by iteration until the wake geometry con- 
verges. The influence coefficients include, 1) the near wake, 2) the intermediate 
wake, and the 3) the far wake. The influence coefficients caused by the near wake 
which are not contained in the code, TFAR1, are computed everywhere in the near 
field of the blade. These influence coefficients are added to those associated with 
the intermediate and far wakes. By multiplying the final strength of bound circula- 
tion to the added influence coefficients, the induced velocities are computed every- 
where in the near field of the rotor blade. The linear compressibility correction 
is made by including the lift-curve slope change in the circulation-downwash rela- 
tion. The bound circulation and blade lift for a given blade geometry are calcu- 
lated by iteration until the wake geometry, and hence the bound circulation, no 
longer changes during iteration. The wake geometry is initially generated by apply- 
ing the momentum theorem, assuming constancy in the flow field, and by integrating 
the downwash velocities from it. The downwash caused by the wake vortices which are 
not contained in the computational domain of the TFAR1 code is obtained on the 
blade, that is, in the near field, from the calculated bound circulation and wake 
geometry . 


RESULTS AND DISCUSSIONS 


The correct position of the wake is very important to obtain correct pressure 
distributions. Existing finite-difference codes usually capture the near wake which 
is presented by the jump in potential across a flat continuation surface traced by 
the blade. The roll-up and nonlinear interaction of wake vortices are not included 
in TFAR1. The roll-up and wake geometry are computed in the free-wake lifting-line 
code, OPLIN, and the velocity field near the blade caused these wake vortices is 
computed and added to the velocity field of the finite-difference code as shown in 
flow chart (fig. 2). The rotor blade geometry used in the computation is shown in 
figure 3. 

Figure 4 shows the pressure distributions from two difference calculations and 
experimental results at spanwise positions r/R = 0.68, 0.80, 0.89, 0.96, tip 
Mach = 0.877 and pitch = 8°. The present results of the coupling of the two pro- 
grams overestimates the pressure, compared to experiments of reference 6 while the 
finite-difference modeling of reference 1 underestimates the pressure. The theoret- 
ical drag coefficient, Cd, is the induced drag due to downwash. All three sets of 
results show the shock formation near the tip, but the calculations capture the 
shocks behind the experimental which may be partially due to disregarding the bound- 
ary layer and partially due to the incorrect downwash field. Figure .5 shows the 
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wake tip-vortex positions from the experiment and the result calculated from the 
OPLIN code. The calculated tip-vortex value is slightly less than the experimental 
value. Figure 6 shows the pressure distributions obtained by the present calcula- 
tion using the coupling of TFAR1 and OPLIN and the experiments of reference 6 at 
five different spanwise locations and at the same rotor blade condition as for 
figure 4. In the inboard portion of the blade, the theory overpredicts the pressure 
distributions and near the blade tip, the predicted shocks are behind the experimen- 
tal locations. 


Figure 7 compares the computed pressure distributions with the measured test 
data at tip Mach = 0.794. Figure 8 compares the pressure distributions at tip 
Mach = 0.612. In both cases, the blade is at a constant pitch of 8°. Both compari- 
sons of calculations and experiments show the same trend shown in figure 6, that the 
different tip Mach number has little effect on the comparison results. Figure 9 
shows the pressure distributions at tip Mach = 0.226 and pitch angle of 12°. 

Except at r/R = 0.5, the agreement is very good, including the suction peak- 
pressure distributions. Figures 10 and 11 show the pressure distributions at tip 
Mach = 0.610, 0.794 for a pitch angle of 12°. The computed shock positions are 
slightly downstream from the experiments, which may be due to the boundary layer 
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vortex beneath the blade tip is necessary to obtain good agreement with experi- 
ments. Figures 12 and 13 show better agreement with experiments with one vortex 
line correction at tip Mach = 0.794, 0.226 and pitch angle of 12°. 
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CONCLUSIONS 


The transonic full -potential flow analysis code, TFAR1, has been modified so 
that the downwash from the code, OPLIN, is added to TFAR1 for hovering performance 
prediction. Also, TFAR1 has been upgraded to compute the flow in the propeller and 
wind turbine cases when the downwash caused by wakes is calculated if the mean 
velocity is normal to the rotordisk. 

The computational efficiency remains about the same as before, but the large 
amount of blade twist required to include the downwash effects is. replaced by the 
additional downwash velocity. The comparison shows that the pressure distribution 
is slightly better-predicted for the higher-loading case of pitch = 12° (figs. 9 
and 10) than for the case of pitch = 8° (figs. 6, 7, and 8); possibly because of the 
smooth wake geometry and downwash field. 
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igure 1.- Geometry of near, intermediate, and far wakes for hovering 
rotors using vortex rings and cylinders to represent the wake, (a) Side view of 
rotor wake model showing intermediate and far wakes formed from vortex spiral for 

2 blades, tip vortex only shown. Blade one Blade two, 

(b) Plane view showing near wake, (c) Formation of intermediate wake. 



Figure 2.- Flow charts for the integration of programs TFAR1 and OPLIN. 
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Figure 5.- Comparison of calculated tip vortex positions with experiments 

of reference 6. 
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- Comparison of present calculated pressure distributions (O) with 
experiments (-o) at tip Mach = 0.877 and pitch angle = 8 °. 
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Figure 12.- Comparison of present calculated pressure distributions (Q-) with 
experiments (o) at tip Mach = 0.794 and pitch angle = 12° (with one vortex line 
correction. 
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